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Using first-principles atomistic simulations, we study the response of atomically-thin layers of 
transition metal dichalcogenides (TMDs) - a new class of two-dimensional inorganic materials with 
unique electronic properties - to electron irradiation. We calculate displacement threshold energies 
for atoms in 21 different compounds and estimate the corresponding electron energies required to 
produce defects. For a representative structure of M0S2, we carry out high-resolution transmission 
electron microscopy experiments and validate our theoretical predictions via observations of vacancy 
formation under exposure to a 80 keV electron beam. We further show that TMDs can be doped 
by filling the vacancies created by the electron beam with impurity atoms. Thereby, our results not 
only shed light on the radiation response of a system with reduced dimensionality, but also suggest 
new ways for engineering the electronic structure of TMDs. 

PACS numbers: 31.15.es, 61.80.Fe, 61.72.Ff, 68.37.Og, 81.05.ue 



Isolation of a single sheet of graphene in 2004 pQ in- 
dicated that strictly two-dimensional (2D) materials can 
exist at finite temperatures. Indeed, inorganic 2D sys- 
tems such as individual hexagonal BN and transition 
metal dichalcogenide (TMD) layers were later on man- 
ufactured by mechanical [2 , 3: and chemical [4[ [5] exfo- 
liation of their layered bulk counterparts, as well as by 
chemical vapor deposition [6j [7] . Recently, TMDs with a 
common structural formula MeX2, where Me stands for 
transition metals (Mo, W, Ti, etc.) and X for chalcogens 
(S, Se, Te), have received considerable attention. These 
2D materials are expected to have electronic properties 
varying from metals to wide-gap semiconductors, similar 
to their bulk counterparts [HJ [9] , and excellent mechan- 
ical characteristics [10]. The monolayer TMD materials 
have already shown a good potential in nanoelectronic 
EJ E] and photonic [1 H31 E] applications. 

Characterization of the fr-BN [T5HT7] and TMD 
[6] [18] samples has extensively been carried out using 
high-resolution transmission electron microscopy (HR- 
TEM). During imaging, however, energetic electrons in 
the TEM can give rise to production of defects due to 
ballistic displacements of atoms from the sample and 
beam-stimulated chemical etching [19 , as studies on h- 
BN membranes also indicate [T5HT71 [20] . 

Contrary to /i-BN, very little is known about the effects 
of electron irradiation on TMDs. So far, atomic defects 
have been observed via HR-TEM in WS2 nanoribbons 
encapsulated inside carbon nanotubes at electron accel- 
eration voltage of 60 kV [21 as well as at the edges of 
M0S2 clusters under 80 kV irradiation [22], while no sig- 
nificant damage or amorphization was reported for M0S2 
sheets at 200 kV [18 - a surprising result taking into 



account the relatively low atomic mass of the S atom. 
Clearly, precise microscopic knowledge of defect produc- 
tion in TMDs under electron irradiation is highly desir- 
able for assessing the effects of the beam on the samples. 
This knowledge would allow designing experimental con- 
ditions required to minimize damage, as well as devel- 
oping beam-mediated post-synthesis doping techniques. 
Moreover, information on the displacement thresholds is 
important in the context of fundamental aspects of the 
interaction of beams of energetic particles with solids, as 
the reduced dimensionality may give rise to an irradia- 
tion response different from that in the bulk counterpart 
of the 2D material [23] . 

Here, by employing first-principles simulations, we 
study the behavior of a representative number of TMDs 
(21 compounds) under electron irradiation, and calcu- 
late the threshold energies for atomic displacements in 
each system, as well as displacement cross sections as 
functions of electron beam energy. In the case of M0S2, 
we also carry out HR-TEM experiments and provide ev- 
idence of electron-irradiation-induced production of va- 
cancies in this material. In addition - inspired by the 
recent advances in introducing impurities in /i-BN mono- 
layers [24j [25J - we discuss irradiation-mediated doping 
of TMD materials. 

For all calculations in this work, we rely on the 
density-functional theory (DFT) with the PBE exchange- 
correlation functional [26] and the projector augmented 
wave formalism as implemented in the simulation pack- 
age VASP [27] [28]. In order to obtain a comprehensive 
picture of the irradiation response of TMDs, we consider 
a large set of layered TMDs: MoX 2 , WX 2 , NbX 2 , TaX 2 , 
PtX 2 , TiX 2 , and VX 2 (where X=S, Se, or Te), which 
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have similar crystal structures. 

We started our study by calculating the displacement 
threshold energy T d (the minimum initial kinetic energy 
of the recoil atom) for sputtering an atom from the ma- 
terial. As in our previous simulations for graphene [29] 
and BN [20] monolayers, an initial velocity was assigned 
to the recoil atom (corresponding to instantaneous mo- 
mentum transfer from the electron to the atom during 
the impact), then DFT molecular dynamics was used to 
model the time evolution of the system. In practice, the 
initial kinetic energy of the recoil atom was increased un- 
til it was high enough for the atom to be displaced from 
its lattice site without an immediate recombination with 
the resulting vacancy. The calculations were carried out 
using a 5 x 5 super cell of a MeX2 monolayer. Test simula- 
tions for larger systems gave essentially the same results. 
The atomic structure of a M0S2 layer and the simulation 
setup are shown in Fig. [TJa) . 

Td required for displacing a chalcogen atom from the 
bottom layer of the sheet [cf. Fig. [IJa)] are presented in 
Fig. [2] In addition to the prototypical MoX 2 and WX 2 , 
we also present results for TiS 2 and TiTe 2 . As evident 
from the figure, Td G [5, 7] eV for all studied compounds. 

We also calculated the vacancy formation energies (Ef) 
for each of the compounds to see how it correlates with 
T d . We defined E f as: 



Ef = E vac - (E hu \ k — Mx)? 



(i) 



where Ekuik and E vac are the energies of the pristine and 
vacancy containing supercells, respectively. The chemi- 
cal potential fix of the chalcogen species is taken as the 
energy of the isolated atom to enable a straightforward 
comparison with the results of dynamical simulations. 
Ef, with and without relaxation of the atomic structure 
of the layer, is also presented in Fig. [2] 

In the non-relaxed case, the energetics is very similar 
for all materials, and - as can readily be noticed - the 
agreement between Td and the non-relaxed Ef is strik- 
ing. This is because during the sputtering of chalcogen 
atoms from the outermost layer, little energy is transmit- 
ted to the surrounding metal atoms due to sufficiently 
fast sputtering event and the rigidity of the structure. 
The energies for the relaxed geometries show a more in- 
triguing behavior. Atomic relaxation for some systems 
evidently gives rise to a considerable drop in Ef, so that 
the similarity to Td is lost. This drop quantifies the de- 
gree of structural relaxation around the vacancy, which 
is minor for MoS 2 , see Fig.[TJb). The analysis of the elec- 
tronic structure revealed an occupied bonding type va- 
cancy state close to valence band maximum and an empty 
anti-bonding type state in the mid gap, which stabilizes 
the structure. For the occupied bonding defect state, the 
electronic charge is localized at the vacancy site, anal- 
ogous to bulk MoS 2 , where Mo atoms donate electrons 
to S atoms. This is true for all of the semiconducting 



materials: MoX 2 , WX 2 , and PtX 2 . The rest of the con- 
sidered materials are metals or semimetals, for which the 
bonding vacancy state may become unoccupied, which is 
reflected in larger relaxation and lower formation ener- 
gies. 




FIG. 1. (Color online) (a) The setup used in the dynamical 
DFT simulations of atom sputtering from TMDs under elec- 
tron irradiation. Initial energy acquired due to the impact of 
an energetic electron was assigned to the recoil atom, then 
DFT molecular dynamics was used to model the evolution of 
the system, (b) M0S2 sheet with an S vacancy. The charge 
densities of the occupied and unoccupied defect states are 
visualized by (blue) transparent and (red) solid isosurfaces, 
respectively. 
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FIG. 2. (Color online) Displacement threshold energies Td 
obtained from DFT molecular dynamics calculations (crosses) 
and formation energies of chalcogen vacancies with non- 
relaxed (dashed lines) and relaxed (solid lines) geometries in 
transition metal dichalcogenides MeX2, X=S, Se, Te. 

Knowing Td, it is possible to estimate the electron 
threshold energy through the relativistic binary collision 
formula and the atom displacement cross section (for rel- 
atively light atoms) by using the McKinley-Feshbach for- 
malism [30]. Ef for Se and Te compounds are smaller 
than in S compounds, but - due to the higher atomic 
mass - their creation through ballistic electron impacts 
requires significantly higher electron energies. In the case 
of MoS 2 , MoSe 2 , and MoTe 2 , T d of 6.9, 6.4 and 5.9 eV 
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correspond to electron energies of about 90, 190, and 270 
keV, as calculated assuming a static lattice. For other 
compounds the required electron energies should be of 
similar magnitude, based on the close values of in 
Fig. [2] An accurate estimation of the displacement cross 
section requires including the effects of lattice vibrations 
on the energy transferred from an electron to a target 
atom [3T] , We calculated the cross sections for vacancy 
production as a function of electron energy for M0S2, 
WS2, and TiS2 beyond the static lattice approximation, 
as shown in Fig. [3] We stress that the production of S 
vacancies for practically all TMDs is within the energies 
commonly used in TEM studies. 

The displacement thresholds for chalcogen atoms in 
the (top) layer facing the beam proved to be consider- 
ably higher than for the bottom chalcogen layer, as the 
displaced atom is "stopped" by the other layers. How- 
ever, after a vacancy is created in the bottom layer, the 
threshold energy for the top S atom in M0S2 to be dis- 
placed and fill the vacancy is about 8.1 eV. This is similar 
in magnitude to the threshold for displacing S atom from 
the bottom layer (6.9 eV), and thus formation of vacancy 
columns should be possible even at 80 kV when lattice 
vibrations are accounted for. for transition metals are 
even higher, since they are bonded to six neighbors and 
similarly stopped by the S layer. For instance, about 
20 eV is required to displace Mo atom from its site in 
the M0S2 lattice, which corresponds to electron energy 
of 560 keV. Naturally, under such conditions the S sublat- 
tice is quickly destroyed. Formation of transition metal 
vacancies is thus considered highly unlikely. 

With regard to possible vacancy agglomeration under 
continuous irradiation, we found that creation of a va- 
cancy does not alter the formation energy in the neigh- 
boring sites in the semiconducting TMDs. Thus, we do 
not expect accelerated formation of large vacancy clus- 
ters. In the same vein, however, it is worth noting that 
chalcogen atoms may also be sputtered fairly easily from 
the edges of nanostructures [2TJ|22]. For example, our 
calculations for a WS2 ribbon show that the chalcogen 
atoms at the edge can have a displacement threshold as 
low as 4.2 eV, as compared to 7.0 eV away from the edge. 

To check our theoretical results on irradiation-induced 
vacancy formation in M0S2, we experimentally studied 
the evolution of a M0S2 sheet under an 80 keV electron 
beam. First, free standing single layer M0S2 samples 
were prepared by mechanical exfoliation of natural M0S2 
bulk crystals, followed by characterization via optical mi- 
croscopy on a Si+90nm Si02 substrate and transfer to 
a perforated TEM support film (Quantifoil), similar to 
graphene samples [32]. The TEM grid was adhered to 
the Si02 surface by evaporating isopropanol on top of it. 
After this, the silica was etched with KOH. Aberration- 
corrected (AC) HRTEM imaging was carried out in an 
image-side Cs corrected FEI TITAN microscope at a pri- 
mary beam energy of 80 keV. The contrast difference 
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FIG. 3. (Color online) Cross section for sputtering a sul- 
fur atom from M0S2, WS2, and TiS2 sheets as calculated 
through the McKinley-Feshbach formalism and the dynami- 
cal values of the displacement thresholds. Dotted lines are 
the data for the static lattice, and solids lines are the results 
of calculations where lattice vibrations are taken into account 
assuming a Maxwell-Boltzmann velocity distribution. The 
cross denotes the experimentally determined cross section for 
MoS2- The inset shows the same data for a larger range of 
electron energies. 

between the Mo and S sublattice is clearly detectable in 
the AC-HRTEM images proving the single layer nature of 
the sheet (for a double layer the contrast would be iden- 
tical as Mo is stacked above S). This is also confirmed in 
diffraction measurements, as successive diffraction spots 
from one {hkl} family show different intensity, whereas 
for bi- and multilayers they are equal [18] . The analyzed 
intensity ratio of the {1100} diffraction spots was found 
to be 1.07. 

During continuous imaging we found an increasing 
number of vacancy sites (exclusively on the S sublattice) 
accompanied by crack formation [see Fig. |4ja)] and lat- 
eral shrinkage of the membrane. Counting the actual 
number of sputtered atoms as in Ref. 31, the cross sec- 
tion for sputtering was found to be 1.8 barn, which is in 
a reasonable agreement with the calculated cross-section 
of 0.8 barn, taking into account that the theoretical es- 
timates are very sensitive to inaccuracies in the parame- 
ters of the model (e.g., and the velocity distribution) 
at energies below the static threshold. 

In Fig. [4^d,e) we present simulated TEM images [33] 
for the single and double vacancies, respectively, based 
on atomic structures (Fig. |4jb,c)) obtained from the 
DFT calculations. Similar defects are observed in the 
experimental TEM images [Fig. gf)]. Different defects 
can clearly be distinguished by analyzing the (Michelson) 
contrast relative to the contrast of the Mo atoms in the 
pristine area. We find that the experimental (simulated) 
ratios are 0.9 (0.9) for a sulfur column, 0.5 (0.4) for the 
single- and 0.2 (0.2) for the double vacancy. 
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FIG. 4. (Color online) AC-HRTEM images of single- layer 
M0S2. Atoms are dark and white spots correspond to the 
holes in the hexagonal structure. During continuous 80 keV 
electron irradiation a crack formed in the membrane (a). The 
edges are terminated with Mo atoms, as can be determined 
from the contrast in the inset. Also an increasing number of S 
vancancies is observed. Most of them are single vacancies but 
also double vacancies, where the top and bottom S atoms are 
removed, can be found. Structure models and corresponding 
HRTEM image simulations for both vacancy types are shown 
in (b, d) for the single- and (c, e) for the double vacancy. 
Experimental examples are shown in (f). To enhance the 
visibility of the defect, a Gaussian filter (0.7 A FWHM) was 
applied to (d, e) and (f). Scale bars are 5 nm (a) and 5 A (f). 



FIG. 5. (Color online) (a) Substitution energies of the impu- 
rity atoms in M0S2 layer calculated for three different chem- 
ical potentials: isolated atoms (blue), diatomic molecules 
(red), and molecules with hydrogen (green). The solid phase 
reference is also given for carbon and silicon (black), (b) Re- 
spective local density of states of the impurity atoms, the 
spin- up and spin-down components. The gray areas denote 
the band edges of the MoS 2 layer, (c-e) Series of AC-HRTEM 
images demonstrating vacancy filling. The red arrow high- 
lights an initial S vacancy that picks up an atom between 
(d) and (e), and the green arrow indicates an S atom that is 
sputtered away between (c) and (d), forming a single vacancy. 



Having shown that vacancies can be created in TMDs 
under electron irradiation, we move on to study whether 
they could be consecutively filled with other atomic 
species deliberately introduced into the TEM chamber. 
We calculate the formation energy of substitutional de- 
fects in M0S2 and consider donors F, CI, Br, and I; ac- 
ceptors, N, P, As, and Sb; double acceptors C and Si; 
hydrogen H and H2; and isoelectronic species O, S, Se, 
and Te. The formation energies and the local density of 
states (LDOS) around the substitution site are shown in 
Fig. [5j We list the formation energies with three different 
chemical potentials of the substituted (impurity) species: 
the isolated atoms, diatomic molecules, or molecules with 
hydrogen (CH 4 , SiH 4 , NH 3 , PH 3 , AsH 3 , HF, HC1, and 
HBr) where we set /in = ^^h 2 - 

Due to the high formation energy of the vacancy, all 
substitutions are energetically favored with respect to 
the isolated atom. However, with respect to \i in the 
molecule, C, Si, and N substitution have positive forma- 
tion energies. Obviously, even if the equilibrium energet- 
ics does not favor the formation of the substitutional de- 
fect, the substitution may still be achieved under electron 
beam, because molecules like N2 or hydrocarbons will 
constantly break apart under the electron beam. Thus 
post-synthesis electron-mediated doping may also be re- 
alized in 2D TMDs, similar to BN sheets [11 [25]. The 
LDOS shows that N, P, As, and Sb behave as acceptors, 



whereas F, CI, Br, and I are likely to be donors. C and Si 
have levels in the middle of the gap and the isoelectronic 
species like O, Se or Te do not produce any localized 
states, as expected. 

Filling of the vacancies was also observed in the TEM 
images, as shown in the series of panels in Fig. [5jc— e) . 
Although we could not identify the type of the impurity, 
this example proves that electron-beam-mediated dop- 
ing is possible. Consequently, through control of atomic 
species in the TEM chamber and the choice of the elec- 
tron energy, modification of the physical properties of 
TMDs via electron beam should be attainable. 

To conclude, we calculated atom displacement thresh- 
old energies in a number of TMDs. These energies are 
a measure of the radiation hardness of the material, and 
serve as critical input parameters in the Kinchin-Pease 
and other semiclassical theories of defect production and 
ion stopping [23j [34] . Here we use them to calculate elec- 
tron displacement energies and corresponding sputtering 
cross sections to quantitatively assess the amount of dam- 
age created in 2D TMD materials during a TEM exper- 
iment via knock-on processes. Observations of vacancies 
in our experimental AC-HRTEM images of single M0S2 
sheets validate our theoretical predictions. Finally, we 
observe filling of the vacancies and discuss the prospects 
for electron-beam mediated doping of TMDs. 
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